Echo spectroscopy of bulk Bogoliubov excitations in trapped Bose-Einstein 

condensates 



m ■ 
o ■ 
o : 

m ■ 
m ■ 



> 

oo 
in 

&\ 
O 

m 
o 



c 

o 
o 



13 



E. Gershnabel, N. Katz, R. Ozeri, E. Rowen, J. SteinhauerQ and N. Davidson 
Department of Physics of Complex Systems, 
Weizmann Institute of Science, Rehovot 76100, Israel 

We propose and demonstrate an echo method to reduce the inhomogeneous linewidth of Bogoli- 
ubov excitations, in a harmonically-trapped Bose-Einstein condensate. Our proposal includes the 
transfer of excitations with momentum +q to — q using a double two photon Bragg process, in which 
a substantial reduction of the inhomogeneous broadening is calculated. Furthermore, we predict an 
enhancement in the method's efficiency for low momentum due to many-body effects. The echo can 
also be implemented by using a four photon process, as is demonstrated experimentally. 



Bragg spectroscopy of a trapped Bose-Einstein conden- 
sate (BEC) has recently revealed many of the BEC's in- 
triguing bulk properties such as global coherence 0] , ver- 
ification of the Bogoluibov excitation spectrum, superflu- 
idity and the superfluid critical velocity, and suppression 
of low momentum excitations due to quantum correla- 
tions of the ground state However, whereas the 
resonance frequency and the frequency integral of the 
dynamic structure factor S(k,u>) mainly reflect the bulk 
properties of the BEC, the frequency width of the spec- 
trum is dominated by inhomogeneous broadening mech- 
anisms 0. The inhomogeneous broadening is due to 
Doppler broadening and the inhomogeneous mean 
field energy that is often well-described within a local 
density approximation (LDA). Even when the LDA pic- 
ture is not complete, as in the recent observation of radial 
modes within the condensate, the inhomogeneous broad- 
ening (manifested as the envelope function of the multi- 
mode spectrum) still dominates @0- Overcoming the 
inhomogeneous mechanisms (arising from the particular 
trap shape), opens the possibility of studying the ho- 
mogeneous broadening mechanisms, that can reflect the 
intrinsic decoherence processes of the bulk excitations, 
e.g. elastic collisions with the BEC ||. 

In this letter we propose a novel echo spectroscopic 
method, which reduces the inhomogeneous broadening 
of bulk Bogoliubov excitations in a BEC (in analogy to 
the familiar echo in two level atoms [jj]). This method is 
based upon the transfer of Bogoliubov excitations from 
momentum +q to —q, by a degenerate two-photon Bragg 
transition induced by an optical lattice, with wavcnum- 
ber 2q. Using the Gross-Pitaevskii equation (GPE), we 
calculate the line shape of the echo excitation spectrum 
and show a substantial reduction of the inhomogeneous 
broadening. We also show a surprising quantum enhance- 
ment of the efficiency of this method for low momentum, 
in contrast with the familiar suppression due to structure 
factor considerations Q. Finally, we demonstrate the 
echo concept experimentally using an alternative scheme 
based on a four-photon Bragg transition. Such a scheme 
allows us to create +q excitations by a two-photon Bragg 
process Q Q and then transfer them to — q excitations via 
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FIG. 1: The echo mechanism, (a) By application of the echo 
standing wave (Sui « 0), +q excitations are transferred to — q 
excitations. The echo process can be achieved either by using 
a degenerate two-photon process (solid arrow) or a degenerate 
four-photon process (dashed arrow), (b)-(d) Cross sections of 
the wave function in longitudinal momentum space during the 
various echo steps, for a momentum q r = 8.06/im -1 , calcu- 
lated with the GPE. The radial direction has been integrated, 
(b) After a 1.3 msec adiabatic Bragg pulse, (c) after the Bragg 
pulse and 4.2 msec delay, (d) As in (c), but with an adiabatic 
4.2 msec two-photon echo pulse instead of the delay. 



a degenerate four-photon Bragg process, using the same 
laser beams. 

As shown in Fig. ^d, a short two-photon Bragg pulse 
with wavenumber +q initially excites the condensate uni- 
formly 10], along the z axis of the cylindrically symmet- 
ric condensate. We then employ our echo scheme, shown 
schematically in^,. Specifically, the initial Bragg pulse 
is followed by a second, much longer Bragg pulse with 
wavenumber 2q. For the frequency difference Slo ~ be- 
tween the two Bragg beams, the echo resonance condition 
is fulfilled for a transition from the positive momentum 
q to the negative momentum —q, as indicated in Fig. 
by the solid arrow and Fig. ^i. The various dispersion 
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curves in Fig. are a schematic representation of the 
different local Bogoliubov dispersion relations due to the 
inhomogeneous local mean-field. They can also repre- 
sent different radial modes The echo line-shape is 
expected to be free of inhomogeneous effects, as long as 
there are no transitions between the various curves, which 
belong to different positions within the condensate (es- 
sentially within the LDA) or different radial modes. Of 
course, the excitations do move along z during the pro- 
cess and therefore every curve is coupled with its vicinity, 
reflecting the fact that q is not a good quantum number 
in an inhomogeneous system of finite longitudinal size. 
This combination of inhomogeneous mean field and fi- 
nite longitudinal size yields a residual small broadening 
of the spectrum. 

To verify the echo concept and calculate this residual 
broadening, we use a simulation of the GPE and 
calculate Bragg and echo processes for our experimental 
system jjj. Our experimental system contains N = 10° 
atoms, confined in a harmonic trap with axial and radial 
angular frequencies w z — 2ir x 26.5 Hz and w r = 2ir x 226 
Hz, respectively, chemical potential fj,/h= 2.06 kHz, and 
LDA average healing length £ = 0.23/im. All excitations 
are axial, maintaining cylindrical symmetry 0|. 

Fig. nj) shows the calculated momentum distribu- 
tion after an initial 1.3 msec Bragg pulse, where the 
q = ground state component and +q excitations are 
seen. The Bragg excitation fraction is given as the in- 
tegral of the wavefunction in momentum space, divided 
by hq. Fig. ^ shows the calculated momentum distribu- 
tion after an additional 4.2 msec of free evolution in the 
harmonic trap. The observed momentum components 
smaller than +q arc the fraction of the excitations which 
have left the condensate bulk during the free evolution 
and were consequently slowed by the external harmonic 
trap. Fig. |2Ji shows the calculated momentum distribu- 
tion after an initial 1.3 msec Bragg pulse and 4.2 msec 
echo pulse, where the appearance of a — q momentum 
population is seen. The echo excitation fraction in Fig. 

is obtained by subtracting its total momentum from 
the total momentum of Fig. ^j, and then dividing by 
2hq. To suppress non-linear effects, we adiabatically in- 
crease and decrease the intensity of our pulses during the 
echo steps. 

Using these GPE simulations, we first compare the 
spectrum of the Bragg and echo excitations in the high- 
momentum, free-particle regime shown in Fig. for 
q r = 8.06fim~ 1 (q£ i = 1.85). Fig. shows the cal- 
culated Bragg spectrum (using a 4.2 msec rectangular 
pulse). The multi-peak structure, corresponding to the 
recently observed radial modes is evident in the spec- 
trum. The calculated FWHM of the spectrum is 1.26 
kHz. The dashed line in Fig. 0} is the LDA line-shape 
0. Although lacking the multi peak structure, the LDA 
line-shape has a very similar width to the exact GPE 
line-shape |7J. Fourier broadening for a 4.2 msec rectan- 
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FIG. 2: Bragg and echo spectra for q = 8.06/im _1 (free par- 
ticle regime), (a) Excitation spectrum of the Bragg process 
calculated with an exact GPE simulation (solid line) and in 
the LDA (dashed line), (b) Excitation spectrum of the echo 
scheme from a GPE simulation. 



gular Bragg pulse, Doppler broadening and the collisional 
broadening (not included in the GPE calculations) are 
210.9Hz, 135.5Hz and 35.1Hz (FWHM) respectively^ all 
much smaller than the inhomogeneous broadening |l2j. 

Fig. |2t> shows the calculated echo spectrum. A narrow 
peak is seen close to Slo = 0, as expected. Comparison 
with Fig. [2K indicates that the echo has indeed reduced 
the inhomogeneous linewidth. Thus, the echo pulse does 
not mix the various radial modes. The echo FWHM is 
0.59 kHz, 2.1 times smaller then the Bragg FWHM of 
Fig. The Doppler broadening, however, is still not 
resolved. The observed ^0.5kHz negative shift of the 
resonance from Slo — is probably due to the decrease in 
momentum of the outgoing excitations during the echo 
process (as explained above). 

We repeat these calculations for low-momentum ex- 
citations in the phonon regime with q — 3.10/Ltm _1 
(q£ = 0.71). The resulting echo spectrum width (Fig. 
Eb) has a FWHM of 0.39 kHz, 1.7 times smaller than the 
FWHM of the Bragg spectrum shown in Fig. Ot- The 
shift of the echo resonance from zero, seen in Fig. |3t>, 
is smaller than the high-momentum case of Fig. |2t>, as 
expected. 

We verify numerically that the width of the echo spec- 
trum increases for shorter and stronger echo pulses, indi- 
cating an increase of Fourier and power broadenings, and 
decreases for longer condensates (having the same /z). 
Thus, having chosen sufficiently weak and long pulses, 
the echo spectra of figures |21 and El are finite-size broad- 
ened. 

Next, we calculate the echo transition rate in the ap- 
proximation of an infinite, uniform Bogoliubov gas, and 
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FIG. 3: Same GPE and LDA calculations as in Fig. for 
the Bragg (a) and echo (b) spectra, only with q = 3.10^tm~ 
(phonon regime). We observe narrowing of the echo spectrum. 
In a addition we note a smaller shift of the echo resonance, 
with respect to fig. 

in analogy to the calculation of the Bragg rate Q. The 
initial interaction Hamiltonian between the light and 
BEC is, 

H = C y CjalckamSl+n-k-m (1) 

klrnn 

Where C is the coupling constant, c (&) are the pho- 
tonic annihilation (creation) operators and a (<r ) are the 
atomic annihilation (creation) operators. The first tran- 
sition we analyze is that from +q to —q, i.e. from the 
initial state \i) = \n q ,n^ q ; N q , N^ q ) to the final state 
I/) = \n q +l,n- q — l;N q —l,N- q +l) where n q represents 
the number of photons in the +q direction and N q rep- 
resents the number of +q Bogoliubov excitations. Since 
the excitations are axial, we refer only to the magnitudes 
of the momenta. In order to evaluate the transition ma- 
trix element, we must consider the matrix element of the 
Hamiltonian between |i) and |/), and transform the 
atomic operators into Bogoliubov excitation operators. 
We find 

(f\H\i)+ = G^/n^^+l^W^N^ + l{u 2 q + v 2 q } 

(2) 

where u q and v q are the Bogoliubov amplitudes [l3J and 
the "+" notation represents the transition from +q to 
—q. The second transition to be explored is the re- 
verse process (transfer from — q to +q), whose matrix 
element is denoted by (f\H |z)_. In the Fermi golden 
rule approximation the transition rate is given by Rate — 
f (\(f\H^) + \ 2 - \(f\H'\i)-\ 2 )6(E q - E. q ). Assuming 
the initial condition N q ^> 1, i.e. (N q + 1) w N q , ne- 
glecting the n q N- q term, and assuming a classical laser 



FIG. 4: The echo response R (solid line). The dashed line 
shows the structure factor S q . Note the enhancement of R 
and the suppression of S q for small q. 



field (i.e. n q ^> N_ q ), the rate can be approximated by 
Rate » WC\ 2 N q n. q n q {u\ + v 2 } 2 6(E q - £L 9 ). Hence, 
there is no Bosonic amplification due to the least pop- 
ulated mode (N- q in our case). Considering the term 
{u 2 + v 2 } 2 and using the normalization u 2 — v 2 = 1 and 
the value of the structure factor S q = (u q — v q ) 2 we find 
that the rate per excitation with momentum q is propor- 

/ i+s 2 \ 2 

tional to the response R, defined by R = ( 2g q J . R has 

the low q asymptotic behavior of (2S q )~ 2 , thus for low q 
(where S q is small) we calculate a large enhancement in 
R. 

The response is seen in Fig. 0] (solid line) to be greatly 
enhanced for small q, in contrast with the Bragg process, 
which is suppressed by the structure factor (dashed line) . 
We emphasize that the validity of these results is for the 
infinite homogeneous gas. We do not observe such an 
enhancement in the GPE calculations for a trapped BEC. 

To demonstrate the echo spectroscopy experimentally, 
we use the apparatus described in in which a nearly 
pure (> 90%) BEC of 1 x 10 5 87 Rb atoms in the 
\F, rrtf) — |2, 2) ground state, is formed in a cylindrically 
symmetric magnetic trap. A two-photon Bragg transi- 
tion and the subsequent four-photon echo transition are 
both induced by the same counter-propagating (along z) 
laser beam pair, locked to a Fabri-Perot cavity line, de- 
tuned 44GHz below the 5S 1/2 ,F = 2 — ► 5P 3/2 ,F' = 3 
transition. 

After exciting approximately 35% of the condensate by 
means of a short 0.1 msec rectangular Bragg pulse with 
Slu = 2n x 16 kHz, we apply a 1 msec echo pulse with an 
envelope shaped as sin(7r x 10 3 t) (0 < t < 1 msec), using 
the same beams. We vary 5uj to make a spectroscopic 
measurement of the echo response around Slu — 0, us- 
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FIG. 5: (a) GPE predictions. The solid curve and the lower 
axis show the four-photon echo spectrum for the smoothly 
varying 1 msec pulse used in the experiment. The dashed 
curve and the upper axis correspond to the Bragg spectrum 
for the same pulse, (b) Experimental echo spectrum (squares) 
and a Gaussian fit (solid line). 



ing the four-photon process to transfer excitations from 
+q to — q (Fig. ^l, dashed arrow). From the absorption 
images after 38 msec of time-of- flight we find the ratio be- 
tween the — q population N^ q and the total (not collided) 
amount of atoms N q + Nbec + N- q . The measured echo 
spectrum is shown in Fig. \Sjp together with a Gaussian 
fit to the data, centered close to zero frequency, and with 
a FWHM of 1.21 kHz. The experimental echo duration is 
limited by collisions, which are not taken into account in 
our GPE simulation and cause a significant degradation 
in our signal for echo times longer than 1 msec. Slosh- 
ing of the condensate may also cause a broadening of 
the measured spectrum width. However, we still observe 
a sub-Fourier broadening of the echo spectrum, which 
is also predicted by a GPE simulation (for four-photon 
echo), as shown by the solid line in Fig. |SJi (FWHM 0.86 
kHz) . The dashed line in Fig. corresponds to a Bragg 
spectrum of a 1 msec pulse (with the same envelope as 
the echo), and is Fourier limited (FWHM 1.63 kHz). The 
sub- Fourier width of the echo spectrum is possible due to 
the non-linearity of the four-photon process. Thus, even 
with relatively short pulses, a narrowing of the lineshape 
is achieved, as compared to standard Bragg spectroscopy 
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In conclusion, we implement a spectroscopic echo 
method in momentum space in order to reduce the mea- 
sured inhomogeneous linewidth of Bogoliubov excitations 
in a trapped BEC. We note that by varying the shape of 
the trap (either in length or in functional form) the colli- 



sional and finite-size broadenings may be reduced further 
than reported here. We also predict an enhancement of 
the echo rate at low momentum emerging from construc- 
tive interference of the amplitudes of various quantum 
paths for this process. This is in contrast to the usual 
suppression of low momentum processes in quantum de- 
generate Boson systems. 
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